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Summary 



This research studied the temperature behavior in the wellbore of the horizontal water 
injector, using the reservoir simulation with the multi-segment well model It examined 
the applicability of the analytical interpretation methods to estimate the flow profile from 
the temperature data The result of this study indicates that the temperature data during 
early injection time and during shut-in period can be used for estimation of the injection 
profile along the wellbore. 

Temperature during early injection time (within few hours immediate after beginning 
of the injection) can be used for estimation of the injection profile with simple analytical 
method, provided enough temperature contrast between injected fluid and the fluid inside 
the wellbore. At a later time of injection, when the temperature profile stabilizes, the 
temperature data gives qualitative information of the injection profile. 

Analytical method for interpreting the temperature profile during shut-in period gives 
detail flow profile along the wellbore. The shut-in analysis for the well with short 
injection history gives better estimation of the injection profile compared to the one with 
long injection history. 

With the enormous data that permanently installed Distributed Temperature Sensing 
(DTS) can produce from a horizontal injection well, identification of the method of 
analysis, the applicability of the analysis method, and the type of data for analysis are 
very important This study gives some insight into and the answer to these questions. 
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Introduction 



The temperature logging has been used long time in the past for profiling the injection 
rate in the vertical well. The analysis techniques for temperature data were based on the 
temperature profile in the wellbore during shut-in period following the injection. During 
the injection, the temperature in the wellbore at the injection interval depends mostly on 
the convection heat rate (the heat carried by fluid). The thickness of the injection interval 
is short, thus the temperature is essentially constant and equal the temperature of the fluid 
at the top of the injection interval (Kg. 1). 




Depth 



Fig. 1 - Temperature in the vertical well during injection 

After shut-in, heat is transferred from undisturbed region of the reservoir to the fluid in 
the wellbore mostly by conduction (Fig, 2). This process leads to the increasing of the 
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temperature in the wellbore and is called the temperature recovery from injection. The 
rate and the magnitude of the temperature recovery depend on the total heat exchange 
between formation and fluid at each location along the wellbore, which is determined by 
the injection rate into reservoir at that interval. Temperature at the interval, which 
accepted more injected water will rise slower than other intervals. Analyzing the rate of 
the temperature recovery in the wellbore allows inferring the injection profile of the 
well. 1 * 5 




Fig. 2 - Temperature profile in the vertical well during shut-in 

Existing methods of analysis of injection profile target the vertical well, where the 
injection interval is usually small and the effect of finite wellbore volume can be 
neglected Further more, the displacement process of the reservoir fluid can be depicted 
by radial flow model with reasonable accuracy. For horizontal well, the displacement 
process occurred along long wellbore. The flow regime starts with radial (ifk^k^ky) or 
elliptical (if h is much smaller than fe, ky) and changes to linear flow in the direction 
normal to wellbore. At the later time of injection, the effect of heat conduction from the 
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cap and base formation can mask the heat conduction in the horizontal direction. The 
finite volume of the wellbore in the horizontal section also contributed to possible 
deviation from existing method of injection profile analysis. All these factors combined 
can make the behavior and the interpretation of temperature in horizontal injector far 
different from that in the vertical well. 

The growing number of the horizontal wells and the development of the new technology 
to accurately measure the temperature in the well without well intervention (DTS) 
revitalize the interest in the temperature logging. Because of the lack of the theory and 
experience in interpreting the temperature in horizontal well, there is a need for a solid 
understanding of the temperature behavior in the horizontal well and to establish the 
procedure for using and interpreting the temperature data from horizontal injection well. 
This study addressed these problems. 

The objective of this study is to: 

- Gain insight of the temperature behavior in the horizontal water injector, and identify 
the difference in temperature behavior between vertical and horizontal well. 

- Investigate the use of the DTS data for the injection profiling in horizontal water 
injector, specifically: the method for analysis temperature data and the type of 
temperature data for analysis. 
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Approach 



I choose to approach this problem by using the thermal reservoir simulation with the 
option of calculating the temperature at any point along the length of the well. 6 * 7 The 
reservoir simulation, which conserves the mass and energy balances of the well-reservoir 
system, allows investigating the effect of different factors on the temperature profile of 
the well. These factors such as the anisotropy and 3D dimension of the injection process 
otherwise could not be captured by the analytical solution. Following stages have been 
performed to achieve the goal of the study: 

- building the simulation model to study the temperature behavior in the well and the 
factors that affect the temperature profile. 

- reviewing the existing interpretation techniques and assessing the applicability of 
these methods for horizontal well using simulated and real temperature profiles. 

- using the simulation model to evaluate the applicability of the interpretation 
technique in analyzing the temperature data from the actual well. 
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Key Results 



The key results from this study are: 

- The temperature data during early injection can be used for estimation of the injection 
profile. This required high frequency temperature along the wellbore within few 
hours of injection. 

- At a later time of injection, when the temperature stabilized, the temperature profile 
provides qualitative information of the injection profile in the well. 

- The shut-in data gives more detail infonnation of the injection profile than the data 
during injection. The new analytical method, which assumes two distinctive 
temperature regions in the reservoir at the shut-in gives good estimation of flow 
profile compared to the simulated flow result The shut-in analysis for the well with 
short injection history gives better estimation of the injection profile compared to the 
one with long injection history. 

- Simulation is a good tool for justifying the estimated flow profile and for 
incorporating other variables (such as pressure distribution, formation permeability 
distribution) into the analysis. In using simulation to investigate the temperature 
behavior in well, it is important to choose the right grid scheme and grid size around 
the wellbore. 
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Simulation Study 



This section briefly described the simulation model and the temperature behavior of the 
well during injection and shut-in. The model simulates the single horizontal water 
injector in a bounded reservoir. Special attention has been paid to the selection of the 
right block size around the wellbore. The specific features of the temperature profile in 
homogeneous and heterogeneous reservoir are discussed. 

Simulation Model 

Fluid Model. The simulator used in this study is compositional Eclipse-300. The fluid in 
the reservoir is modeled by two components, one is CI (40%), the other is C7+ (60%). 
The mixture of these two components under the reservoir condition used in this study 
results into an undersaturated oil with GOR of 300 scf?STB and bubble point pressure of 
2500 psi. 

Reservoir model. The reservoir is of 1 10 ft thick with the aerial extend of 7000 ft in X 
direction, 4100 ft in Y direction. It is filled with oil and irreducible water saturation of 
0.2. Porosity of the reservoir is 02 and kept unchanged throughout the course of the 
study. The thermal capacity of the rock is constant The thermal conductivity and the 
permeability of the reservoir are subject to change during the course of the study. Table 1 
shows the basic parameters used in this simulation study. 

TABLE 1 - BASIC RESERVOIR PARAMETERS 



Reservoir thickness, ft 


110 


Permeability, md 


3 


Porosity 


0.2 


Initial bubble point pressure, psi 


2500 


Irreducible water saturation 


0.2 


Geothennal gradient, °F/ft 


0.02 


Temperature at top of the formation, °F 


210 


Rock heat capacity, Bta/ti?-°T 


35 


Rock thermal conductivity, Btu/ft-D-°F 


24 



Well model. The well is horizontally completed at the center of the reservoir in the X 
direction and is modeled as a number of segments with the ID of 0.45 ft The pressure 
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drop in each segment is calculated using homogeneous flow model The heat transfer 
coefficient between wellbore segment and reservoir is calculated based on the formation 
thermal conductivity between the wellbore and the well block. In doing so, we assumed 
that the fluid inside the well segment have the same temperature as calculated at the 
sandface. The well parameters and its thermal heat resistance coefficient are given in 
Table 2. 



TABLE 2 - WELL PARAMETERS 



Well ID, ft 


0.45 


Wellbore roughness, ft 


0.001 


Well length (open hole), ft 


5000 


Wellbore heat resistance coefficient, l/(Btu/ft-D-°F) 


0.008 



Temperature Behavior 

Grid selection. The temperature profile in the wellbore is affected by the fluid 
convection and the heat conduction between the wellbore and reservoir. Since the 
temperature behavior of the well depends on the temperature distribution around the 
wellbore, refined grid block scheme must be used The first task of simulation study is to 
build a model that can capture the true temperature behavior in the wellbore with a 
reasonable computational cost One approach to compromise the accuracy and the 
computational cost is to choose the model with an asymptotic temperature behavior. In 
this method, we refine the grid dimension of the model, starting from the coarse grid and 
study the temperature behavior during injection and shut-in. We chose the model, starting 
from which the more refinement of the grid scheme does not notably change the 
temperature behavior of the well. 

The first model has uniformed grid block size in all direction: 100 ft in X direction, 
100 ft in Y direction and 10 ft in Z direction. At time t=0, reservoir is at the initial 
geothermal temperature. Water is injected with constant rate for 3 days and shut in for 3 
days (Fig. 3). The temperature along wellbore is reported from simulation at each time 
step during injection (Fig. 4) and shut-in (Fig. 5). 
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Fig. 3 - Well performance 
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Fig. 4 - Temperature profile at different time during injection. 
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Fig. 5 - Temperature profile at different time during shut-in. 



As I refined the grid size in Y and Z direction, the calculated temperature at well 
segment significantly changes, especially during shut-in. Figs. 6-7 show the temperature 
of the well segment at the middle of the well at different time during injection and shut-in 
for different grid refinement schemes. We can see that the temperature during shut-in 
changes very much as we refine the grid size around the welibore. The temperature 
profile stabilizes as the size of the well block is 2 ft in Y, Z direction and exponentially 
increases outward (cases 4-6 on Fig. 7). The model, starting from which the temperature 
does not change with further refinement has 19 grid blocks in Y direction and 13 layers 
as shown in Table 3. 

TABLE 3 - GRID BLOCK DIMENSION 



1024, 512, 256, 128, 64, 32, 16, 8, 4, 2, 4, 8, 16, 32, 64, 128, 256, 512, 1024 



Y,ft 



Z,ft 



22,12,8, 5,4,3,2,3,4,5, 8, 12, 22 
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Fig, 6 - Temperature at the middle of the well during injection for different models. 
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Fig. 7 - Temperature at the middle of the well during shut-in for different models. 
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To make sure that this model results a good temperature behavior with the change of 
rate and injection time, another model with excessive refined grid scheme in the Y and Z 
directions has been built This model consists of 79 grid block in Y direction and 2 1 grid 
blocks in Z direction. This model and the model shown in Table 4 were run at two 
different injection rates (one is 5000 STB/D, another is 10000 STB/D). The temperature 
at the middle of the well shows very similar behavior with time (Fig. 9). 



180 
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Refined Model, Q2 



Q1= 5000 STB/D 
Q2 = 10000 STB/D 



Fig. 9 - Temperature behavior at the middle of the well: selected model and refined 
model with different rates. 



For homogeneous reservoir, the grid block size in X direction does not affect the 
simulated temperature in the wellbore after the temperature stabilized. But it does affect 
the calculated temperature along the wellbore at the early time of injection. For 
heterogeneous case, the grid block size in X direction is 25 ft - arbitrarily chosen to 
capture the heterogeneity along the wellbore. 
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The simulation model shown in Table 3 is used to gain insight of the temperature profile 
of the well. The well performance is similar to the one shown in Fig. 3. We will first look 
at the temperature behavior in the homogeneous reservoir, and then in the heterogeneous 
reservoir with different injection profile. In the real life, the temperature of the injected 
water at the heel of the well is a function of time. In order to accurately represent the real 
temperature behavior during early time of injection, we specified injection temperature 
according to the observed temperature from a well with similar configuration and 
injection rate (Fig. 10). This temperature of the injected fluid at the heel of the well 
during early injection time can be estimated using Ramey approach. 



250 




0 0.5 1 1.5 2 2.5 

Time, day 

Fig. 1 0 - Injection temperature at the heel of the well. 

Homogeneous reservoir. At very early time, we can observe the advancing of the 
injection front toward the toe of the well (Fig. 4). The higher the injection rate or the 
smaller the wellbore diameter is, the fester this front is advancing. With the injection rate 
of 5000 STB/D in homogeneous reservoir for the well having ID of 0.45 ft and length of 
5000 ft, the simulation shows that the injection front will advance to 70% of the well 
length for 1 hr. At a later time, the temperature profile is stabilized with little change 
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along the wellbore. Time for the temperature to be stabilized depends on the heat thermal 
conductivity of the reservoir and the injection rate. For the injection rate used in this 
study, it takes 2 days for the temperature to be stabilized in a homogeneous reservoir. 

During shut-in period, the temperature recovers in the wellbore. The rate and magnitude 
of the temperature recovery depend on total injection volume, the thermal conductivity of 
the well (Fig. 9). For the well in homogeneous reservoir, the rate of temperature recovery 
is constant for the whole length of the well (Fig. 5). 

Heterogeneous Reservoir. Similarly to the temperature behavior in the homogeneous 
reservoir, the temperature profile can be considered as differently behaves in 3 periods. 
During early injection, the velocity of injection front depends on the distribution of 
injection flow along the wellbore. Fig. 11, 12 show the temperature profiles for two 
different distributions of the permeability along die wellbore. The temperature chest 
moves faster at the interval with low permeability and slower in the interval with high 
permeability. 




Distance along Wellbore, ft 



Fig. 1 1 - Temperature profiles at early time of injection: permeability increases along 
wellbore. 
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Distance along Wellbore, ft 

Fig. 12 - Temperature profiles at early time of injection: permeability decreases along 
wellbore. 

At a later time, temperature profile only bears the information of the injection profile for 
the end of the well. If the injection rate at the end of the well is high, the temperature 
profile is flat (Fig. 13). If the injection rate at the end of the well is low, the temperature 
profile is curving upward (Fig. 14). 
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Distance Along Horizontal Section, ft 



Fig. 14- Stabilized temperature profile for the case with low injection rate at the tip of 
the well. 
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For all cases of the permeability distribution, the shut-in temperature profiles correlate 
with the injection profile along wellbore (Figs. 15-17). The rate and the magnitude of the 
temperature recovery are higher in the interval with small injection rate than in the 
interval with high injection rate. 




Distance Along Horizontal Section, ft 



Fig. 1 5 - Temperature profile during shut-in, low permeability at the middle. 
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Distance Along Horizontal Section, ft 



Fig. 1 6 - Temperature profile during shut-in, low permeability at the tip. 




Distance Along Horizontal Section, ft 



Fig. 17 - Temperature profile during shut-in, low permeability at the heel. 
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Effect of Rate. The high injection rate case differs from low injection rate case in 2 
areas: it accelerates the rate of cooling in the wellbore due to larger volume of cool water 
being injected into reservoir; and it creates higher pressure drop in the wellbore or change 
the flow profile along wellbore. Three cases have been run to study the effect of the rate 
and pressure drop along the well bore on temperature profile: one with injection rate of 
5000 STB/D, one case with injection rate of 1 5000 STB/D; and the third case with 
injection rate of 15000 STB/D and reduced wellbore diameter (from 0.45 ft to 0.35 ft). 
The total pressure drop from the heel to the toe of the well is 4, 50 and 200 psi 
respectively. Fig. 18 compares the injection temperature profile for the first two cases. 

» 
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Fig. 1 8 - Temperature profile during early injection: different injection rate. 

The only difference in the temperature profile is mat the temperature profile for high 
injection rate quickly stabilizes along the wellbore. The injection profile is almost 
.constant for the case with small pressure drop in wellbore. In the case of high pressure 
drop in the wellbore, the injection profile slightly decreases along the wellbore 
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(maximum difference of up to 1 0% - Fig. 19). During shut-in, the magnitude of the 
temperature recovery is very small in the case of higher injection rate and it does not 
capture the change of the injection profile along the wellbore (Fig. 20). It is worth noting 
that the shape of injection profile along the wellbore is different, depending on the 
position of the well in the X direction. The injection profile would have the U shape if the 
well locates far from the edge of the reservoir. The result implied by the temperature 
behavior as shown in Fig. 20 is that temperature profile may not capture the small • 
variation of the injection profile along the wellbore. 
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Fig. 20 - Temperature profile during shut-in: injection rate Q = 1 5000 STB/D 



For the model with high pressure drop in the wellbore, a notable fluid redistribution may 
occur in wellbore after shut-in (Fig. 21). This phenomenon may affect the temperature 
profile along the wellbore. As we can see, the region with higher injection rate will have 
fluid coming back into wellbore during shut-in. 
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Fig. 21 - Flow rate to/into reservoir after the shut-in (negative value means flow from 
reservoir, positive means flow into reservoir) 



Conclusions 

- The selection of the grid block size of the model has great effect on the 
temperature behavior in the wellbore. The grid block size is as small as 2 ft 
around the wellbore is required for simulating the temperature profile of the well. 

• Compare to the vertical well, the temperature profile during early injection time is 
controlled by the injection profile and injection rate. . 

- Temperature profile at late injection time may give only qualitative information of 
the injection profile. 

Temperature profile during shut-in is controlled by injection profile, thermal 
conductivity of the reservoir, injection rate and injection time. 
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Interpretation Methods of Temperature Profile 



This section describes the interpretation technique and their application. Three 
interpretation techniques are considered: one for interpreting the temperature profile 
during the early injection period based on the movement of the water in the wellbore; one 
for shut-in temperature data based on the source and sink method for temperature build 
up in the wellbore; one for shut-in temperature data based on the assumption of two 
temperature region in the reservoir at die shut-in. 

Analysis of Injection Data 

In the previous section, the temperature behavior has been studied for different flow 
profile. There is no interpretation technique for interpretation of the temperature during 
injection for vertical well Because of the short wellbore length at the injection interval, 
temperature in wellbore quickly stabilized to the temperature of the injected water. For 
horizontal well, however, if the temperature along the wellbore before the injection is 
known, the temperature profile can be used to estimate the flow profile, providing that 
enough contrast between wellbore temperature and the temperature of injection fluid 
exists. Assuming that we can determine the speed at which water moves inside the 
wellbore from the temperature profile, we can determine the flow profile at any interval 
along the wellbore as: 




where the subscript i denotes the location of interest, 7 is running from the heel of the well 
to the location of interest As we can see from die equation, the accuracy of injection 
profile depends on the accuracy of the total injection rate, the wellbore ID, and the 
position of the injection front in the wellbore. 
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Analysis of Shut-In Data 

Nowak method. Nowak 1 proposed the method for determining the injection profile using 
the approximation of the source and sink method for temperature build-up in the vertical 
wellbore. This results to a linear relationship between injection intake rate and the 
temperature recovery along the wellbore at any time / after shut-in: 

^ = A to r77 < 2 > 

where AT) is the temperature difference between the reservoir temperature and the 
temperature at time /after shut-in; Q is the total heat transfer in the reservoir, which can 
be assumed proportional to the product of injection rate and injection time qfa . Subject 
to a constant injection rate, the Eq.2 can be expressed as: 

flW=^A». (3) 

This equation is used in analysis the injection profile along the wellbore as followed: 

1- The shut-in temperature curve in the region of non-permeable strata above the 
injection interval is extrapolated through the injection interval for determining the 
reservoir temperature. 

2- The area between the actual curve and the extrapolated curve is determined. 

3- For each unit depth, the fraction of the total injection rate is the same as the fraction 
of the area between two curves over the total area. 
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Depth 

V 

Fig. 22 - The Nowak method for injection flow profiling. 

The Nowak method does not account for the heat transfer in the vertical direction. Thus, 
the transition zone between injection interval and impermeable strata can be erroneously 
interpreted (Fig. 22). Cocanower et al? proposed a method for better determining the 
boundary between formations. This method used temperature profile at early after shut-in 
to calculate the temperature decay coefficient and used this coefficient to calculate the 
shut-in temperature at a later time. However, this method is affected by the fluid 
redistribution in the wellbore at early time after shut-in as discussed in previous section 
and hence will not considered here. 

New alternative method. Assuming that the temperature distribution in the reservoir at the 
shut-in has the shape of two distinctive regions, one with average reservoir temperature, 
one with the temperature of the wellbore at the shut-in (Fig. 23), the temperature 
behavior at the wellbore can be expressed as followed: 2 

r.-^JfrS/.g)* (4, 
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Fig. 23 - Assumption for the temperature distribution in the reservoir: two temperature 
regions at the shut-in. 




Dimensionless Time, at/Rc 



Fig. 24 - Solution to the temperature behavior at the well in the reservoir as shown in 
Fig. 23. 
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The solution of this equation is shown in Fig. 24 in term of the dimensionless 
temperature and dimensionless time. Assuming that the temperature front is the same as 
the saturation front, this solution can be used to relate the total injection water into each 
interval and the radius of the temperature front as followed: 

(5) 

Or, in term of total injection rate: 

Qbn = ^EMw , (6) 

;=i 

This equation is used to estimate the flow profile along the wellbore with t D! determined 
from Fig. 24. The different between this method and Nowak method is that it related 
injection rate to the temperature in the exponential, instead of a linear, relationship. 

The procedure for estimating the injection profile is as followed: 

1- Determine the initial temperature of the formation T r . For horizontal well, it is the 
initial temperature of the reservoir. 

2- The dimensionless temperature for each interval is determined from the temperature 
log and the reservoir temperature: Ta=(Tr-T s )/(T r T^). 

3- For each unit depth, determine the dimensionless time t D from Fig. 24. 

4- This dimensionless time for each interval along wellbore is used in Eq. for estimating 
the correction coefficient A. The flow rate into reservoir for this interval is: A*t Di . 



Both aforementioned methods for shut-in temperature data determine the injection profile 
based on the assumption of known reservoir temperature and, small radius of the 
injection front in the reservoir, and no heat transfer in the X direction. For the well with 
long injection time or unknown well history, when these assumptions are no longer valid, 
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the result of analysis may be suffered. In the next section, these methods of estimating the 
injection profile will be examined on the actual well data. The estimated flow profile 
from each method will be compared to the simulated flow profile to see which one 
performs best 
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Application to Safah-200 Water Injection Well 



This section examined the temperature profile of the actual water injection well. The 
numerical simulation model is built to reproduce the measured temperature data* The 
injection profile from simulation is used to compare with the injection profile, estimated 
by using the analytical methods. The model is also used to assess the applicability of 
analytical method at different times of the well life. 



WeU Performance and Temperature Behavior 

This well is completed in a low permeability micritic limestone reservoir. The reservoir 
fluid is light (42 deg. API), low viscosity (0.4 cp). Average porosity is 22%. The 
permeability of the reservoir is less than 5 md The reservoir thickness is of 100 ft -130 
ft. 8 Kg. 25 shows the schematic of the completion of the well. 9 hi the horizontal section 
of 5400 ft, a stinger of 2*7/8 is used to convey DTS along the wellbore. The well is under 
injection with the rate of 4500 STB/D for 14 days, following by 1 day of shut-in and 3 hr 
of re-injection (Fig. 26). The temperature data is available for the shut-in period and for 3 
hr of injection that followed. 
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, Fig. 25 - Well completion: blue line represent the fiber optic line. 
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Fig. 26 - Well injection history. 

After 14 days of injection, the temperature in the well shows that the first half of the well 
length has the almost constant temperature. Temperature at the second half is gradually 
increases to the temperature of the reservoir. The temperature profile is very similar to 
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that for heterogeneous reservoir, where permeability is gradually decreases toward the 
well (Fig. 12). During shut-in, temperature increases at different rate along the wellbore 
(Fig, 27). We will applied the analytical method to estimate the injection profile from 
temperature profile during shut-in. One important observation is that during shut-in, the 
temperature at above casing shoe is increases at a higher rate, compared to the part of 
wellbore below casing shoe. After one day of shut-in, the temperature difference at 
casing shoe is about 40 F. When the injection is resumed, this portion of hotter water is 
pushed along the wellbore (Fig. 28). The water temperature is decreasing partly due to 

the water going into formation, partly because of the heat loss while moving toward the 

* i 

toe. The movement of this hot slug of water along the wellbore can be used to estimate 
the injection profile along the wellbore in the same manner as described in previous 
section for temperature during early injection. 
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Fig. 28 - Temperature profile during three hours of re-injection. 

In the next section, we present the result of the injection profile estimation from three 
analytical methods as described in the previous section for temperature data during shut- 
in and re-injection. 

Shut-in Analysis 

Nowak method Following the procedures mentioned in previous section, the flow profile 
is estimated from temperature data during shut-in. Fig. 30 shows the injection profile 
estimated from the temperature at the end of the shut-in period by using Nowak methods. 
Note that for constant initial reservoir temperature, the injection profile is just the curve 
of (T r Ti\ scaled down so that the area under the curve is the total injection rate. 
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Fig. 30 - Injection profile from Nowak method. 



New method. Assuming the initial reservoir temperature of 226 F, the dimensionless 
temperature at two different shut-in times are shown in Fig. 31. We can see that around 
the tip of the well, the dimensionless temperature becomes unstable because the 
temperature at the shut-in is closed to the reservoir temperature. Any drift or noise in 
temperature may cause the to become unreasonable large or small. 
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Fig. 31 - Dimensionless time at two different shut-in times. 



We can overcome this problem by combining the magnitude of the temperature recovery 
into account of Multiplying the T D with the Nowak temperature (T r T s ) will dampen 
the effect of noise in the region with very small temperature difference. The injection 
profiles using these normalized T D are shown in Fig. 32. 
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Fig, 32 - Injection profile using new method 
Injection Analysis 

Assuming that the position of the slope of the temperature curve in the front of the 
transition zone represents the injection front (Fig. 33), the velocity of the water inside the 
wellbore is shown in Fig. 34. The noise in the velocity data indicates that the position of 
the injection front in the wellbore is sensitive to the noise in the temperature data. From 
Eq. 1, it is clear that the accuracy in determination of injection profile is affect by the m- 
situ total injection rate, the accuracy of front position from temperature data, and 
wellbore diameter. Wellbore ID can be assumed to be constant in well with liner 
completion. The fluctuation of the total injection rate and the accuracy of the injection 
front from temperature data are the major factors that affect the result Because of the 
noise in temperature data, this interpretation method is difficult to be done automatically. 
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Fig, 33 - Example of injection data analysis. 
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Fig. 34 - Velocity of the water in the wellbore, determining from temperature profile 
during early re-injection 
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Assuming constant well ID and constant injection rate, the fluid velocity can be 
represented by a best-fit line (or lines) as shown in Fig. 34. The injection profile can be 
estimated from that fluid velocity as is shown in Fig. 35. 




Fig. 35 - Injection profile for the data in Fig. 34 from temperature during injection. 



The estimated injection profiles from all three methods are shown in Fig, 36. As we can 
see, all methods show similar feature of the injection profile: the first half of the well 
contributes most of the injection rate of the well, while the second half of the well 
contribute small portion of the injection rate. The interpretation method for shut-in data 
gives more detailed feature of injection profile than the method for injection data does. 
For shut-in temperature data, the new method gives the injection profile with higher 
contrast between two parts of the well, compared to the Nowak method. 
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Fig. 36 - Comparison of the injection profile for different method. 



Since we do not have the injection log to justify the accuracy of the analytical method, 
we used numerical simulation to assess the applicability of interpretation methods. We 
used the simulation model described in the previous section to simulate the temperature 
profile of the well. The temperature profile is matched with the measured temperature 
profile by manually adjusting the thermal conductivity of the reservoir and the 
permeability along wellbore. In assessing the applicability of the analytical method of 
flow profiling, we assumed that the best method would have the injection profile similar 
to the simulated injection profile. 



Simulation Study of Safah-200 

The data for the simulation is in Table 4. The reservoir dimension in the direction along 
wellbore is of 21 8 blocks. The well is completed in cell 3 to 216 in X direction, each cell 
has the length of 26 ft Two first blocks and two last blocks' in X direction have the size 
of 50 ft and 100 ft each to represent die effect at the tip of the well. The well is subject to 
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an injection rate of 4500 STB/D for 14 days, and then is shut-in for 1 day. The 
temperature at the heel of the well during early time of injection is the function of time 
and is taken from the similar well Safeh-179 and is shown in Fig. 10. The initial value of 
the thermal conductivity of the reservoir is assumed to be 60 Btu/ft-D-F. The analytical 
estimate of the injection profile from alternative method is used to determine the 
distribution of the reservoir permeability along the wellbore. We assumed the linear 
relationship between flow rate and permeability: krqi. The permeability estimated from 
flow rate for each interval is applied in all grid block in Y and Z directions for that 
interval. With these input parameters, the simulated temperature profile at the t=14 days 
of injection and t=l day of shut-in is shown on Fig, 37, together with measured 
temperature profile. We can see that the temperature profile from simulation is in good 
agreement with the measured temperature in the second half of the wellbore but too high 
for the first half of the well length. Reducing the injection rate in the first half of the well 
would result into smaller temperature recovery for this half of the well, but lower 
temperature in second half of the well during injection. Thus, the good match between 
simulated temperature and measured temperature can only be achieved both by reducing 
the thermal conductivity of the reservoir and by adjusting the permeability profile. 

After several adjustments, the model gives good match with the measured data (Fig, 38). 
The final model has the reservoir thermal conductivity of 30 Btu/ft-D-F. 
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Fig. 37 - Simulated temperature - first run with the injection profile estimated from new 
method and reservoir thermal conductivity of 60 Btu/D-ft-deg.F. 
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Fig. 38 - Simulated temperature - final run with the adjusted permeability along wellbore 
and reservoir thermal conductivity of 30 Btu/D-ft-deg.F. 
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The flow profile from the final simulation model is plot in Rigs. 39-40 for comparison 
with the flow profiles, analytically estimated from temperature data. Injection profile 
estimated with the alternative method for shut-in temperature data is in very good 
agreement with the profile from simulation. The Nowak method underestimated the 
injection flow rate in the first half of the well. It is interesting to note that both methods 
for shut-in data overestimate the low flow rate and underestimated the high flow rate. 
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Fig. 39 - Injection profiles from different methods. 
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Fig. 40 - Comparison of the injection profile between analytical analysis and simulated 
result 



Effect of Injection History on Shut-in Analysis 

While the analysis method for injection temperature data does not involve the complex 
heat transfer phenomenon and can be applied at any time during the life of the well, the 
shut-in temperature data depends on the injection time and the initial reservoir 
temperature. For the well with long injection time or unknown injection history, the 
applicability of this analysis method is not certain. We tried to take a further step toward 
the study of how this method behaves with time. For this purpose, the simulation model 
obtained by matching the temperature with observed temperature is run for 275 days. To 
simulate the possible-effect of the heat transfer from top. and base formations, 4 additional 
layers with zero permeability are introduced to the model Two layers on the top have the 
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thickness of 50 ft and 100 ft. Two layers at the bottom have thickness of 50 ft and 100 ft. 
The heat capacity of those four layers is set at 350 Btu/ft3-F. 

During 9 months of injection, five shut-in periods are taken place at 1 day, 5 days, 1 5 
days, 100 days, and 270 days (Fig. 41). 
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Fig. 41 - Well rate history for investigation the temperature behavior with time. 

The plot of the temperature profile at the time before each shut-in period is shown in 
Fig. 42. The temperature profile during injection shows that the temperature front 
continues advancing toward the toe of the well (which is consistent with the actual data 
from other wells). However, the injection profde along the wellbore shows very little 
change with time (Fig. 43). It suggests that the change of the temperature profile during 
the injection after stabilization does not indicate the change in flow profile, but can be the 
effect of heat conduction process. This observation, however, does not exclude the 
possibility of the change of the injection profile in the real field condition, when the 
operation of other wells changes the pressure distribution in the reservoir. It only suggests 
that the change of the temperature during injection at the late time can be contributed 
mainly by heat conduction. 
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Fig. 42 - Temperature profile at different time of injection. 
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Fig. 43 - Simulated injection profile at different injection times. 
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The temperature profile at the end of each shut-in period is used for estimation of 
injection profile with the new analysis method. The result of analysis for each shut-in is 
shown in Fig. 44. The comparison with the actual flow profile is shown in Fig. 45. The 
estimated flow profile from shut-in data at different shut-in period shows that the 
accuracy of the analytical technique for injection profile deteriorates with time. This can 
be explained that with the long injection time, two major assumptions of the analytical 
method are being violated. Those assumptions are: 1) no heat transfer in X direction, and 
2) sharp boundary of two temperature regions in the reservoir at the shut-in time. Thus, 
for better estimation of the injection profile at the later time of injection, the method that 
incorporates the history of the well and the heat transfer in X direction is needed. 
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Fig. 44 - Analytically estimated injection profile from simulated temperature for different 
shut-in periods (shown in Fig. 41) 
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Fig. 45 - Comparison between estimated and true injection flow rate at time. 



Conclusion 

The early injection data following short shut-in period can be used to estimate the flow 
profile. This method of flow profiling analysis can be used at any time in the well life and 
only required enough temperature contrast of the fluid in the wellbore before injecting. 
This temperature contrast can be achieved by shutting the well for 1-2 days and letting 
the part of the well located above the injection formation to warm up. However, in the 
case of gradually increase of the permeability starting from the heel of the well, the sharp 
temperature contrast may not be observed. 

The accuracy of the analysis method for shut-in temperature data is decreasing with 
injection time. For the well in this study, the correlation coefficient between estimated 
and true injection profile is slightly decreases with the time. Compare to the analysis of 
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injection data, shut-in data give more detailed flow profile estimation. Table 4 presents 
the comparison between the injection and shut-in analysis of the temperature data. 

Finally, the simulation is a good tool forinvestigating the temperature behavior of the 
horizontal injection well. It also can be used to incorporate other reservoir parameters 
into analysis and for predicting temperature behavior at the later time of injection. 



TABLE 4 - Comparison between injection and shut-in temperature analysis 



Injection 


Shut-in | 


General. Does not depend on injection 
history. 


Depends on injection history. 


Give gross feature of the injection profile 


Give detailed feature of the injection 
profile 


Does not depend on the reservoir properties 


Depend on the near wellbore 
heterogeneity 


Need hi^b frequency injection data- . 
synchronize with total injection rate 


No need for high frequency data. 


Difficult to be done automatically 


Can be done automatically 



For the better estimation of the injection profile in horizontal water injector, more work 
should be dpne to improve the accuracy and reliability of the analytical method, 
specifically in the following areas: 

- determination of the injection front in wellbore from temperature data during the 
early injection. 

i 

incorporation of the well history into analysis of shut-in temperature data. 
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Appendixes 



Shut-in Analysis 

The transient temperature response , of a cool region with the radius of Rc in a hot 
formation can be expressed as: 

a-. 

The solution for r=r w has the form as shown in Fig. 24 and can be approximate as: 
ln(T D )= -1.950071 E" 6 *x 6 +8.955864 E" 5 *x 5 -1.550909E S *x 4 +0.0131941 *x 3 - 

0.05875433 *x 2 - 0.8758939 *x-1.571235 

where x = ln(<w/i? c 2 ), A-2 

Assuming that if the To for some injection interval is known, the value of atlRl can be 

found: T D =f In -p-j, therefore: ln-~ = f'{T D ), or: i? c 2 * ate'''™ 

Then we can expressed the injection rate at this interval as a function of Td as follows: 

e ^{i-sj atfu-sj Mi-sJ 

,Or: 

= pi = gfatie^) = W A-3 



where /' {T D ) is the solution for ln-^- at To as shown in Fig. 24; A can be estimated as 

J!, 



c 



^ .-ryl-Qt A-4 

Since the temperature distribution in the actual reservoir is not the same as assumed, A 
can be used as a regression coefficient for all intervals with the constraint of total 
injection rate: 

Q^A^e-''^, A-5 

i 

The interpretation procedure with the account of the magnitude of temperature recovery 
is following: 

1- Assume reservoir temperature T r and calculate the maximum of (T r - T^) 

2- Calculate the dimensionless temperature T m = r - ' sl r 

3- Solving for $ (To) according to Eq.A-2 

4- Estimate A, Eq. A-4 

5- Calculate Eq. A-3 

6- Adjust A so that total injection rate is match, Eq. A-5 
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Injection Analysis 



Consider the horizontal well with constant injection rate Q. The well has length of L. The 
fluid is incompressible* The wellbore diameter is a function of the distance along 
wellbore and is denoted as D x . The intake rate along the wellbore is assumed to be 
constant with time and is denoted as q(x). Suppose that at time /, the displacement front is 
at distance jc from the heel of the well At time f , the position of the front will be 
x + dbc. The movement of the displacement front is caused by the difference between 
total volume of injected fluid entering the well Qdt and the total volume of fluid going 



X 



into formation in interval from heel to that front dt $q{C)dC 

o 

On 

dV = Qdt-dt)q{C)iC 9 (B-l) 

o 

On the other hand: 

dV = ^dx, (B-2) 

4 

Them 

^•dx = Qdt-dt)q(C)iC, (B-3) 

Since the speed of the displacement front is v c = ^ , we have: 

at 



v,=^rfe-|^j, (B^) 

Eq. B-4 suggests that if we can measure the speed of the displacement front (by 
measuring temperature along wellbore at every time during early injection), we can 
construct the intake profile along wellbore by discretizing the integral term in Eq. B-4. 

^(Xi-x«)=(t,-*^^ (B- 5 ) 

on 

*-i£d tt - s t*-$*>-'4 • *« 

where subscript / denotes die wellbore interval, at which the intake rate being determined. 
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Temperatur Analysis for Horizontal Injection W Us 

Introduction 



Horizontal wells have many advantages over conventional vertical wells, and have 
been extremely popular in the oil industry. The rapid increase in the applications of 
horizontal well technology brought an impetuous development of the procedures to 
evaluate the performances of horizontal wells. 

Many papers have focused on the pressure analysis of horizontal wells. However, 
there are some limitations to the pressure analysis. One limitation is that distributive 
pressure sensing is not practical now. Another limitation is that the conventional transient 
pressure analysis assumes the wellbore is infinite conductive, i.e., the flux along the 
wellbore is uniform. This assumption is not reasonable for long horizontal wells. Pressure 
drop due to friction and non-uniform distribution of formation properties, such as skin 
factor, permeability and porosity, can result in non-uniform flux distribution. So, flux 
distribution is a key factor for horizontal well performance evaluation. 

The DTS, distributed temperature sensing, has been successfully used for different 
wells; The distributive temperature profiles along horizontal wellbores are easily 
measured. The temperature profiles contain useful information of flux distribution. 

The main objective of this paper is to establish a practical model to predict the 
temperature distribution along the wellbore when the injection flux is specified, or to 
estimate the injection flux distribution with measured temperature profile. 



Assumption 

(1) The pressure and temperature gradient in the direction parallel to the axis of the 
horizontal well (jc) are much smaller than those in the directions (y and z) 
perpendicular to the axis. So, the mass and heat transference along x-direction can 
be neglected. 

(2) The viscosity and density of oil and water are constants. So, the pressure and 
water saturation distribution are independent of temperature. 

(3) Isotropic formation. 

(4) Incompressible of injected fluid. 

(5) Open hole completion. 



Wellbore Model 

(1) Wellbore flow rate distribution equation 

Let Q^(x 9 t) denote the injected fluid volume flow rate along the well bore (m* lh) t 
and q^ix^t) denote the fluid volume rate injected into the formation for unit length of 
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the wellbore (m 2 lh).The mass conservation of injected fluid in the wellbore gives (see 
FigJ): 

^=-*„(*>0 (I) 
Integrating equation (1) from the heel to toe yields: 

L 
0 

0 

where fi w (r) is the total injection rate at the heel of horizontal well. 

t 




Fig. 1 Wellbore flow rate distribution 

(2) Temperature distribution equation 

The total heat stored in the tiny element dx is 

=c t AJ w {x,t)dx (3) 

Let q^(xj)=27ir w n^^ denote the heat flow rate from formation into the unit 

length wellbore, T w {x,t) denote the temperature profile along the wellbore. Energy 
conversation equation for element dx is: 

*■«-€..-#- < 4 > 

where and q m:a denote the heat rate flowing into and out of the tiny wellbore 
element dx . a . is composed of two terms: the heat carried by the fluid flowing into the 

^* win * 

element through the wellbore cross-section area at x, c,Q M {x>t)T w (x t i), and the heat 
flowing from formation to the wellbore element through the wellbore surface due to heat 
conduction, q„,(x,t)dx. q^ is also composed of two terms: the heat carried by the 
fluid flowing out of the element through the wellbore cross-section area at x+dx, 
c,Q„,(x+dx,t)T w (x+<ix,t), and the heat carried by the fluid flowing out of the element 
through the wellbore surface, c,g M( (x,0r w (x,/)dx. 

9,*, =c,e irf (x > r)r B ,(x,0+9iw(*.0dx 
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= c,Q« (x + dx, t)T w (x +6x,t)+ c, qv , (x, t)T w (x,t)dx 
Substituting , q m and equation (3) into (4) gives: 

Substituting equation (1) into the above equation yields- 

where c, - heat capacity of injected fluid (J/(mV*)) ; 
4, - flowing area in the wellbore (m 2 ) . 

ij - thermal conductivity of the fonnation(j/(m." KM)), 

Near Wellbore Heat Transportation - Water Injection 

The near wellbore flow regime can be regarded as steady-state radial flow. Let us 
consider a toy radial element between r and r+dr. The total heat stored in this element 
for water injection is composed of three terms: the heat stored in the water phase 
Q„=2m^r^ w ( x ,r,t)cJT(x t r,t), the heat stored in the oil phase 
Q 0 =2wir4{l- Sw (x,r,t)]p 0 T( x ,r,t), and the heat stored in the rock 
a=2^dr(l-^) c ,J(j C ,r,0,i.e., 

Qr«=Q w +Q 0 +Q, 
= 2^drr( JC ,r,/)[ Sw (j: > r,/Mc w -c 0 )+c^+(l-^) Cr ] (6) 
where T(x,r,t) is the temperature distribution in the reservoir, 0 is porosity, c is the 
heat capacity of oil (J/(m\\K)) , and s w (x,r,t) is water saturation. 

dT 

Let q T (x,r,t) = 2mir— denote the heat inward radial flow rate, q„(x,r,t) denote 

the water volume outward radial flow rate, and q 0 (x,r,t) denote the oil volume outward 
radial flow rate. The energy conservation equation for this radial element is- 

'-^- = 9rin-q mtt (7) 
where q m and 9ml denote the heat rate flowing into and out of the toy radial 
element dr . ?rt) is composed of two terms: the heat carried by the ofl and water flowing 
into the element through the inner surface at r, kf.fcr.i)+c A fcr.O]lfer.0.aiid 
the heat flowing into the toy element through the outer surface at r+dr due to heat 
conduction, q T (x,r+dr,t) . is also composed of two terms: the heat carried by the 
oil and water flowing out of the radial element through the outer surface at r+dr 
[c v qjx,r+dr,t)+c 0 q 0 ( x> r+iir,t)]r(x,r+dr,t), and the heat flowing out of the toy 
element through the inner surface at r due to heat conduction, q T (x,r,t). 
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In* = [c w q w (x, r,t)+c 0 q o (jc, r, /)JT(jc, r, i) + q T (x, r + dr, /) 

?« a i=k?w(^'*+<^.0+c^ ( ,(x,r+dr,/)]r(x,r+dr,0+ ?r (r,/) 

37* 

4f r (x,r,/)=2aJ7r-- 
dr 

Substituting the above equations and equation (6) into equation (7) yields: 

[s„(x,r,tMc w -c o )+c 0 t+a-4)c^+Mc„-c o )T(x,rj)^ 

at dt 



1 

2^r 



(8) 



Let V ( x » 0 denote to* water 60111 at tune * » and choose r £ (x, /) . Then we have: 

s w (x,r,t) = l- Sor , q w (x,r > t) = q wl (x,t), 9o (x,r,/) = 0 
Near the wellbore, the heat flux can be also considered as a constant, i.e., 

And equation (8) becomes: 

b-v^-o+^+a-^J?-- 5 ^ 5 ^ (9) 

<w 2^ or 

At the wellbore ( r = r w ), we have: 

fc-'^-o+^+o-^--^f2) do) 

Combining equation (10) with equation (5) and choosing c, = c w yield: 
Let 

«, w =0-^fl-^l+A+(l-<*A, a 2 =^f- 



LA. 

Equation (1 1) becomes: 



I ,w2 ^(x,/)J 



^+Za w (x,,)|k=:0 (12) 

Near Wellbore Heat Transportation - Oil Injection 

For oil injection, the water phase flow rate q w (x,r,t) = 0 and water saturation 
s„(x,r,f) = s wl . The total heat stored in the radial element is composed of three terms: the 
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heat stored in the water phase Q w = 2»tJr^ w c w 7'(jc,r,0 , the heat stored in the oil phase 
Q 0 =2nrdr$l-Srt]p 0 T(x i r,t) ,' and the heat stored in the rock 
a=2^dr(l-f>)c r T(x,r,/),i.e, 

a. 

= 2nrirT(x,r,t)[s^(c w -c o )+cj+(l-0)c r ] 

Similarly, we have: 

9*=c o ? o (x,y)r(x,r,0+? r (*,r + dr,0 
1™ =c 0 q 0 {x,r + dr,t)T(x,r + dr,t) + q r (r,t) 
The energy conservation equation for this radial element is: 
Substituting the above equations into equation (7) yields: 

In the steady state flow regime near the well bore, we have: 

And equation (13) becomes: 

[s«*c v -cJ + c> + (l-^)c r ]^=M ? M^ (14) 

ot 7Jor or 

At the wellbore (r = r w ), we have: 

k**. -o+ «-#+a-«c,]5k« -M*M( K) ( i5) 

Combining equation (15) with equation (5) and choosing c, = c 0 yields: 

^ «,.=^/l-^)+A+(l-^. a 2 =^i 
• «„<x,/)=. ^ 



Equation (16) becomes: 

^• + la a (x,r)^=0 (17) 

Boundary Condition and Initial Condition 

The boundary condition is specified with the injection temperature at the heel 
(x=0): 

^(0,0=7^(0 (18) 
The temperature at the heel T^if) can be determined with the wellbore heat 
transmission model put forward by H. J. Ramey. 
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The initial condition is: 
where T R is the reservoir temperature. 



(19) 



Wellbore Temperature Prediction - Forward Problem 

When the reservoir properties (porosity 0 and permeability it), fluid properties 
(density p 0 andp w , viscosity // 0 and// tf , relative permeability k m mdk m ), wellbore 
geometry (wellbore diameter r w , fluid flow area 4,, roughness e and the length of 
perforated section L) are specified, then the injection rate distribution q^ixj) can be 
determined with analytical model (e.x., model established by TUPREP) or numerical 
model (e.x., ECLIPSE100). And thus, with the properly defined boundary condition (18) 
and initial equation (19), we can predict the wellbore temperature profile T w (x,t) by 
solving equation (12) for water injection or equation (17) for oil injection. 

Let q 0 = QiS^A de note the average injection flux at time / = / 0 , and 



a w =- 



9o 



a.. 



I Wo) 



C *>?0> 



Equations (12) and (17) can be rewritten as 

m 2 

Note the unit of q^fat) and q 0 is 7 — 9 so the unit of a w and a n is 

hour hour 

Let 5 = ~L ' * D ~ °J for water "J^ 011 or t o = f °r oil injection. They denote the 



(20) 
(21) 



1 



dimensionless variables. The dimensionless forms of equations (20) and (21) are: 



Or, both equations can be rewritten as: 



dt r 



(22) 
(23) 

(24) 



ke* £c('d) denote the characteristic curve along which the temperature keeps 



i.e., 
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«-.-l»*iH£-» (25) 

Comparing equation (24) with (25) yields: 

(26) 

Equation (26) is the characteristic equation with respect to the partial differential 
equation (24). Equation (26) defines a group of curves, characteristic curves. We can 
prove that all characteristic curves do not intersect with each other. If one characteristic 
curve crosses the positive g coordinate, then the temperature on this curve is specified by 

the initial condition, Le., equals to the reservoir temperature T R ♦ Otherwise, the curve will 

cross the positive t D coordinate, and the temperature on this curve is specified by the 

boundary condition r w (f^) , where is the intersection of the characteristic curve with 

the time coordinate. 



{ 



> 
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